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Abstract 

We derive an expression for the luminosity distance as a function of redshift for a flat Robertson- 
Walker spacetime perturbed by arbitrary scalar perturbations possibly produced by a modified 
gravity theory with two different scalar perturbation potentials. Measurements of the luminosity 
distance as function of redshift provide a constraint on a combination of the scalar potentials and 
so they can complement weak lensing and other measurements in trying to distinguish among the 
various alternative theories of gravity. 
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I. INTRODUCTION 



General relativity (GR) is ingood agreement with all astrophysical observations of binary 
pulsars and solar system tests jj. These observations provide tight constraints on deviations 
from GR on scales that are smaller or comparable with our solar system. Observations on 
larger scales are less restrictive, so it is possible that gravity is substantially different from 
GR on these scales. In recent years there has been a considerable effort to construct theories 
that modify GR on large scales, partly for the purpose of explaining the current phase of 
accelerated expansion of the Universe without introducing a dark energy component. Among 
;hese theories are scalar tensor theories (see e.g. , f(R) theories and DGP gravity 

Many studies have discussed comparing observables in these theories with present and 
future observations, and in this way constraining and sometimes refuting these alternative 
theories of gravity , lUj). In this paper, we focus on one observable 



- the luminosity distance as function of redshift for a perturbed Robert son- Walker (RW) 
Universe. We calculate this observable for a class of metric theories of gravity including GR. 
Measurements of the luminosity distance as function of redshift for ty pe l a supernovae (SNe) 
provide evidence that the Universe expands at an accelerating rate jlj . . These studies 
would be extended by the planned joint dark energy mission (JDEM). Observations of the 
luminosity distance as function of redshift may be able to constrain the various alternatives 
to GR. 

For any metric theory of gravity the RW metric provides a good description of the Uni- 
verse on large scales. The construction of this metric is based on observations that on large 
scales the Universe is homogeneous and isotropic (see e.g. jl4,[l5|), and also on the assump- 
tion that the Copernican principle holds, namely that we are not in any special location in 
the Universe. Since these general considerations are independent of the theory of gravity 
the construction of the RW metric remains valid for a large class of modified GR theories. 
Assuming the RW metric, one can normally adjust the parameters of a given modified GR 
theory so that it would produce the observed expansion history of the Universe. However, 
small deviations from homogeneity and isotropy which give rise to metric perturbations do 
depend on the particular theory of gravity in use. This makes the perturbed RW metric an 
ideal framework for studying observables that could distinguish among theories of gravity. 

Many previous studies have suggested using cosmic microwave measurements, weak lens- 
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iiiEf measurements and other observations to distinguish amoner the various modified GR 
theories (see e.g. . Recently Bonvin, Durrer, and Gasparini have suggested 



that measurements of the luminosity distance power spectrum induced by cosmological per- 



turbations may be used to determine cosmological parameters [16j. Cooray, Holz, and 
Huterer have showed that two-point angular correlation function of SNe can provide useful 
data to study the foreground of large scale structure 17|. SNe surveys normally have smaller 
data sets then weak lensing surveys so they typically have a larger statistical uncertainty. 
Therefore, with regard to sensitivity to cosmological fiuctuations the SNe surveys are not 
competitive with weak lensing surveys. Nevertheless, they provide an independent measure- 
ment of a different physical quantity (SNe surveys are sensitive to the luminosity of SNe 
while weak lensing surveys are sensitive to the distortion of galaxy images) and so they can 
be used to complement weak lensing measurements. 

Another motivation to study the luminosity distance fiuctuations is that they degrade 



the accuracy of the determination of cosmological parameters from SNe data 
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211]. Estimation of the systematic error produced by cosmological perturbations is therefore 
necessary for SNe luminosity distance data analysis. 

The lack of tight constraints on the theory of gravity on large scales together with the 
sensitivity of future SNe surveys to the foreground of cosmological perturbations motivates 
the calculation of luminosity distance as function of redshift for a perturbed RW Universe 
in modified GR theories. 

In this paper we shall consider a fiat RW metric with linear scalar perturbations. In a 
Newtonian gauge this metric takes the form 

ds'^ = \_-dif{\ + 2^) + (1 - 20)dx • dx] . (1.1) 

Here aij]) denotes the scale factor as function of conformal time, and the potentials and 
satisfy |0(x, ?7)|, \^{y^^r])\ « 1. GR in the absence of anisotropic stresses gives ip = (p. 
The luminosity distance for a perturbed RW Universe have been studied before by several 
authors. Thus Sasaki [22] has studied the luminosity distance as function of redshift for 
a general perturbed spacetime. While Sasaki's analysis is very general, it gives an explicit 
expression for the luminosity distance only for the case of an Einstein-deSitter Universe 
with ip = (f). For the case of = j^an explicit expression for the luminosity distance was 



derived by Pyne and Birkinshaw 



23 ] and was later corrected by Hui and Greene [20|. An 
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equivalent expression was recently derived by Bonvin, Durrer, and Gasparini [16|. In this 
paper we calculate the luminosity distance as function of redshift for the metric (11.1 ). We 
use a different method of calculation than the methods used in Refs. 16, 

y,y,13, and 



our result generalizes the expressions in these references since we allow for possibly different 
scalar potentials (p ^ ip. Our expression fl6.1l) reduces to the corresponding expression in 
Ref. |20|] for the special case of = ■?/' [see Eq. (C21) in this reference]. 

Our complete expression for the luminosity distance Dl a.s function of redshift for the 
metric (11.11) is given by a somewhat cumbersome formula (16.11) . Fortunately, in practice 
most SNe surveys are sensitive only to subhorizon density perturbations. For these surveys, 
we may drop the terms that are subdominant for subhorizon perturbations, and obtain a 
simpler expression for the (subhorizon) luminosity distance Df^^ reading 

Dr\z, n) ^ (;^, - ;^„)(1 + ;.){l + V, ■ n - , ' (1.2) 

^ J Xo Xs Xo "* 

Here the luminosity distance is expressed in terms of the observed redshift z, and the di- 
rection to the source, where n (also denoted as n"") denotes a unit spatial vector from the 
observer to the source. The notation denotes an approximate equality accurate up to 
first-order in the potentials 0, ip, and the peculiar velocities (in the conformal spacetime) 
V. The subscripts s and o refer to the source and the observer, respectively. The conformal 
Hubble rate is denoted Ti = The potentials (p and ip, which by definition are functions 

of the spacetime coordinates, are considered here to be functions of an affine parameter x 
defined on the zeroth order photon null geodesies. The affine parameter x is an implicit 
function of the observed redshift z. This function is determined by the unperturbed RW 
background, and is given by 

- r 1 ^ 

^'~Jo H{z) ^^^"^ 
where H{z) = Ti/a is the Hubble rate, Xo denotes an arbitrary initial value for the affine 
parameter, and we assume that the background spacetime is either expanding or contracting. 

We now briefly discuss Eq. (11. 2p and how it may be used to distinguish among various 
theories of gravity. Notice that this formula has two types of term: terms that are pro- 
portional to the peculiar velocities of the observer and the source, these terms represent 
kinematical Doppler-shift, and a term depending on the Laplacian of the potentials that 
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represents gravitational lensing. This lensing term depends only on the combination ip + (p 
of the two potentials, and so it can not be used to differentiate the GR case ip = (p from the 
more general case (p ^ ip. Fortunately, the velocity terms in Eq. (11.21) provide additional 
information that breaks this degeneracy. The dependence of the peculiar velocities on the 
potentials follows from the equations of motion of the observer and the source. We assume 
that the observer and the source feel no interaction other then gravity, and furthermore we 
assume that the total energy momentum tensor is covariantly conserved. It now follows that 
the observer and the source move along geodesies of the perturbed spacetime, and therefore 
their peculiar velocities satisfy 

V ^ + a,^a"^v + ^ . (1.3) 

Since the peculiar velocities depend on ip, but are independent of (p, they provide information 
that breaks the degeneracy. It is also useful to examine the dependence on the redshift of 
the various terms in Eq. (II. 2p . Terms that depend on peculiar velocities are expected to 
be bounded (the peculiar velocity of a host galaxy is normally of order 500km sec~^) and 
therefore the ratio of a typical peculiar velocity to the Hubble flow is larger at low redshift 
and smaller at high redshift. We might also expect the RMS effect due to lensing to increase 
with redshift. Therefore, the lensing term is expected to be dominant at high redshift, while 
the peculiar velocity terms are expected to become dominant at low redshift. 

Information about the scalar perturbations can be extracted from the data by calculat- 
ing the correlation function [DL{z')DL{z)]^^{DL{zPn')DL{z,n)), where Dl{z) denotes the 
average over angles, and (...) denotes an ensemble average. We anticipate that if we spe- 
cialize to either low or high redshifts, this correlation function would be sensitive to either 
the peculiar velocity terms or the lensing term. This means that by combining information 
from both low and high redshifts one can overcome the above mentioned degeneracy. Such 
measurements may be able to provide information that could distinguish among various 
theories of gravity. 

To calculate the correlation function from theory, one must have some knowledge about 
the underlying theory of gravity. In particular one must know the relation between the 
overdensity and the gravitational potentials (p and tp. Instead of specializing to a particular 
modified theory of gravity it is possible to employ a parameterized framework where the 
relation between the two potentials and the relations among the potentials and the overden- 
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sity are parameterized, such that various modified gravity theories produce different values 

one can use Eq. (11. 2p in 

combination with the power spectrum of density fluctuations to obtain speciflc predictions 
for the power spectrum of the luminosity distance, which can then be compared with ob- 
servations. In this paper we keep our assumptions about the theory of gravity at minimum, 
and we defer the more detailed parameterized analysis to future work. Nevertheless, we can 
still make a rough order of magnitude estimate of the correlation function. It reasonable 
to expect that even in a modifled gravity theory the potentials ip and while not precisely 
equal, should be of the same order of magnitude. This means that while the correlation 



function shou 



d differ from the one in GR, it is likely to be of the same order of magni- 



tude. Ref. [16|] shows that for GR and a CDM Universe the contribution to the correlation 



function form lensing alone can be as large as 10~^ for z = z' = 2 at I 300. Ref. [17| 
estimates that this lensing contribution to the cross-correlation function could be detected 
with a signal to noise ratio of 10 with a survey of 10,000 SNe over lOdeg^ between redshifts 
of 0.1 and 1.7. While the signature of modifled gravity is likely to have a smaller signal to 
noise, it may still be detectable in future SNe surveys. 

This paper is organized as follows. In Sec. [TTl we present the framework for calculating 
luminosity distance in a general spacetime, in Sec. IIIII we specialize to the metric (11.11) and 
simplify the calculation by introducing a transformation to a conformal spacetime, in Sec. 
IIVI we impose initial conditions, and flnally in Sec. IVII we obtain our flnal expression (16. ip . 

II. THE LUMINOSITY DISTANCE IN A GENERAL SPACETIME 



There is a well known general relation 25|, l26| between the observed luminosity distance 



Dl{z) and the observed angular diameter distance Da{z) which reads 

Dl{z) = {1 + z)^Da{z). (2.1) 

This relation is valid in any metric theory provided that the linear momentum of photons is 
conserved, and so in particular it holds for the metric (11.11) . Below we calculate Da{z) which 
is a purely geometrical quantity, and substitute our result back into Eq. (12. ip to obtain the 
luminosity distance Dl{z). 

The angular diameter distance Da{z) is deflned in the following manner. Suppose that 
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an observer views a sizable distant object (e.g. a distant galaxy or a structure of the CMB 
anisotropy) that subtends a small solid angle AQ. Using the geometric optics approximation 
we can describe the electromagnetic radiation using a congruence of null geodesies. We 
consider geodesies that emanate from a vertex at the observer and propagate backwards in 
time towards the source. The angular diameter distance is given by 



Here A denotes the affine parameter along the congruence, and AA{X) denotes the transverse 
cross sectional area of the congruence at a fixed A. 

We pick a representative null geodesic from the thin congruence and denote its worldline 
with x"(A), and denote its tangent null vector field with fc" = where throughout Greek 
indices run from to 3. In general the evolution of a thin null congruence is completely 
described by an expansion parameter 9, a shear tensor cTq,/?, and a rotation tensor Ua/s, where 
the expansion parameter takes the form of 

We assume that the rotation tensor uJajS vanishes at the observer, since this tensor satisfies 
a homogeneous 

everywhere. Therefore, Raychaudhuri's equation for the null congruence reads 

^ = -\(^" - R^pk'^k' -^o.,^""'- (2-4) 
Eqs. (12. 2112. 3112. 4p gives the focusing equation (see e.g. |28l |) 

1 (fDA 

Da dX^ 2 
which we solve to obtain Da- 



--{R^f,k''k'' + a^pa"^), (2.5) 



III. CONFORMAL ANGULAR DIAMETER DISTANCE 

We now specialize the calculation of Da to the metric (II. ip . This metric allows us to 
factor out the dependence of Da on the scale factor, and thereby simplify the calculation. 
To this end we consider the following conformal transformation 

ds'^ = a^ds'^ , 6n = 6n , DA = aDA ,dX = a^dX , k"" = a-^k'' . (3.1) 
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Here quantities with and without tildes denote the conformal space and the real space, 
respectively. Notice that the form of Eq. fl2.5p is invariant under this conformal transfor- 
mation, this means that the focusing equation in the conformal spacetime is obtained by 
adding tildes to all the quantities in Eq. fl2.5p . The focusing equation in the conformal 
spacetime can be transformed into an integral equation of the form 

Da = -\ piRe^ik^'k" + a^p~a^^)DA{X){K - X)d~X + + C^K • (3.2) 

Here Ci and C2 are constants which depend on the initial conditions, and all the quantities 
in the first brackets on the right hand side are evaluated at A. So far we have not used 
a perturbative approximation and Eq. (13. 2p is accurate to all orders in the perturbation 
potentials. Below we calculate Da to the first order in the potentials cj) and ip. 
Expanding the null geodesic worldline a;" (A) in a perturbation series gives 

a;"(A) ^ x°(A) + 5a;"(A) , r (A) ^ fc"(A) + 5A;"(A) , (3.3) 

where an overbar denotes an unperturbed quantity, and a 5 preceding a quantity denotes the 
first-order perturbation to that quantity. In this perturbation scheme the affine parameter 
As characterizes the location of the source on the null geodesic. The same Ag is used both 
for the perturbed quantity and for the unperturbed quantity. Notice, however, that due 
to the perturbation potentials the source is characterized by different coordinates in the 
background spacetime and in the full spacetime ^"(Ag) 7^ a;"(As). Using the above notation 
we also have 

Da^ Da + 5b a , Rap ^ SR^p , a^p ^ 5a^p , Ci^2 ~ C'1,2 + <5C'i,2 • (3.4) 

Notice that at the leading order the conformal spacetime is flat so that the expansions for 
Rap and start at first-order. Using Eq. (13.20 together with Eqs. (I3.3f3.4p we obtain 

Da{\s) = Ci + C2\s, (3.5) 
5bA(\s) = 5Ci + 6C2~Xs ~l [/ SRapk^^DAWiXs - X)d~X . (3.6) 

The flatness of the background spacetime implies that the background null vector k^^ is a 
constant four-vector. By rescaling A we set this vector to be 

/c^ = = (-l,n). 



where n is a unit vector by virtue of the nulhty of k^^. Locating the observer at the origin 
= 0, where a,b = 1, 2, 3, we have 

n{X) =r]o-iX- Xo) , x"(A) = n»(A - A^) . (3.7) 

IV. INITIAL CONDITIONS 

We now supplement equations fl3.5|3.6l) with initial conditions at A = Xo, and thereby 
determine the constants Ci and C2. First we demand that Da{Xo) = 0, which gives 

C*i = ~XoC2 , 6Ci = —X0SC2 ■ 

We therefore have at the leading order Da{Xs) = C'2(As — Xo). Eq. (13.71) implies that C2 = 1. 
With the above initial condition Eqs. (I3.5|3.6p read 

DAiXs) = ~Xs-~Xo (4.1) 
5DaCXs) = 5C2DaCXs) ~l r 5R^pk-¥DACX)iXs - ~X)dX (4.2) 

J Xo 



The constant 5C2 is determined from the value of the derivative at A^. To determine this 

^ dX ° 

constant it is instructive to first calculate where Da is the conformal angular diameter 
distance for a comoving observer, ignoring for the moment the observer's peculiar velocity. 
It is then possible to correct the expression and account for the observer's peculiar velocity. 
In the vicinity of a comoving observer we can Taylor expand the potentials (p aiid by 
treating the distance from the observer as the small parameter. It follows from the metric 
(II. ip that the angular diameter distance takes the form of 



DaCX) = RiX)\/l - 20(A) + OiR'/') . 
where R = (6abx"'x^)^^'^ . Imposing initial conditions Sx^-iXo) = ^^(Xq) = gives 

^(Ao) = l-0o, (4.3) 
ctA 

which gives rise to a constant 6C2 = — 0o for a comoving observer. To correct for the 
observer's peculiar velocity let us consider a transformation to a realistic reference frame 
moving in a velocity Vq with respect to the static observer in the conformal space, where 
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Vo = I Vol = Oi^ip) = 0{(j)). We maintain the notation Da for the angular diameter distance 
in the reahstic frame, which is given by 



Using Eq. fHl2|) together with Eqs. (I4.3II4.4I) we obtain 

J Ao 

Evaluating the expression inside the integral gives 



(5D^(A,) =D^(A,)(vo-n-0,) (4.5) 
/ ' + ^) + 20,^r, - ■ n + - ^,,6)nV] D^(A) (A, - A)dA . 



Using the relation ^ = V0 • n — we integrate by parts and obtain 



(5D^(A,)=5^(A,)(vo-n-0,) (4.6) 
/ ' + V^) - 20,ryr, + (0,a6 - V',a6)nV] 5a(A)(A, - A)rfA - 2 /" ' 0,,(A, + Ao - 2A)rfA . 

J Ao Ao 

V. REDSHIFT 

So far we have calculated the dependence of JDa and on the afiine parameter A. The 
goal of this section is to use these relations to express the luminosity distance as function of 
the observed redshift z. Using Eq. (13.41) and Eq. (13. ip together with Eq. (12.11) we find that 

Dl(A,) ^ (1 + zfa(K) {Da(K) + ^Da(\s)\ . (5.1) 
By definition the observed redshift is given by 

We also introduce a conformal redshift 1 + 5, defined by adding tildes to all the quantities 
in Eq. (15. 2p . Using Eqs. (I5.2|3.ip together with = a~^£t" we find that 

1 + ^ = ^, (5.3) 
a 
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where for convenience we set a(Ao) = 1 at the observer and use the imphcit notation a = 
a{\s). Substituting Eq. (l5.3p into Eq. (15. ip gives 

DLiK) ~ (1 + + z)Da{Xs) + SD^iX,)] . (5.4) 

To obtain the deviation of the luminosity distance as function of redshift due to the metric 
perturbations we need to compare the perturbed luminosity distance and the background 



luminosity distance at the same redshift. Following Hui and Greene 20|] we calculate Dl{Xs + 
6X) the luminosity distance at a shifted affine parameter where the shift 6X is defined by the 
relation 

1 + zCXs + 6X) = 1 + zCXs) , (5.5) 

where 1 + z{f]) = a~^{f]) is the standard redshift in the RW background spacetime. These 
definitions would allow us later to substitute the standard RW relation Xs{z) into Eq. (15.40 . 
Using Eq. (15. 4p and recalling that Da = Xs — Xq we find that 

DlCXs + fx) ^ (1 + z)[DACXs)il + z) + 6X + 6DaCXs)] • (5.6) 

We now calculate z and SX and substitute their expressions into Eq. (15.60 . 

First, we consider z. Using the relation ^ {1 — together with the metric (II. ip 
and the redshift definition (15. 2p applied to the conformal space, we find that 



Xs 



To calculate the quantity (5/c°)^^ we integrate the null geodesic equation. After integration 
by parts we obtain 



5~ (_^ + v-n)J^ - / (0,, + V^,,)dA . (5.7) 

J Ao 

Next, we consider 6X. Combining Eq. (15.31) with Eq. (15.51) we find that 

1 + -zms)] - , (5.8) 

a{r]shift) 

where 

Vshift = vCXs + 5X) + S'niXs) . 

From which we find that 

S~X = Sv-:^, (5.9) 

11 



As 



where we introduced the notation Ti = a^^^a ^ and Si] = [Sx^ was defined in Eq.f l3.3p ]. 
Finally, we calculate 61] by integrating the null geodesic equation which gives 

Sr]C\s) = {(po~i^o)DA+ [ \0,^-^,O(A,-A)rfA + 2 / + 2 / ' ^,^(A,-A)t;A. (5.10) 

To derive this expression we used the initial conditions Sri{\o) = and Sk'^{Xo) = (po + i'o, 
where the last equation is obtained from the requirements k^k^g^y = k^k^rj^y = together 
with the initial condition 5k(Ao) = 



VI. RESULTS 

We now substitute Eqs. fl4.iy4.6p into Eq. ([5S]) and use Eqs. fl5.7y5.9ll510D after some 
integrations by parts and a change of notation, = A, we finally obtain 



dx 

Xo /v,o -J Xo X^ 



Dl{z, n) ^ ixs - Xo)(l + ^){l - + V, ■ n - 2 + f 

Jxo Xs Xo Jx 

~ {Xs-Xo)-Hs (-'P + ''-^)x:- J^^(^^v + tri)dx 

1 f^' , - Xo)(xs - x) , , 1 r%, . I dx-Xo){xs-x) ^ 

V + dx+- / {(f),r,^ + llj,r,rj) dx 

^ Jxo Xs Xo ^ J Xo Xs Xo 



+ 



Xs rx 



-i^,v)dx+ I i<P,v - ^,v) ^° _ ^ dx- ]:i(ps - + (po - i^o)] ■ (6.1) 



Xo Jxo Xs Xo 



This expression gives the luminosity distance as function of the observed redshift and the 
direction to the source for the perturbed RW metric (11. ip . Notice that the last three terms 
in the curly brackets vanish for the case of GR where (p = ip, but may differ from zero for 
modified GR theories. Roughly speaking, the various terms in Eq. (16. ip can be interpreted as 
representing the following physical phenomena. Terms depending on the velocities represent 
kinematic Doppler shift arising from the peculiar velocities of the observer and the source, 
terms depending on the potentials (without derivatives) represent gravitational redshifts, 
and terms depending on a single time derivative of the potentials are analogous to the 
integrated Sachs- Wolfe (ISW) effect. In fact the ISW effect appears in the final term of Eq. 
(15. 7p . It has been previously shown [8] that the ISW effect in modified GR theories produces 



a modification to the low multipoles of the CMB anisotropy power spectra. 

Consider next the first term in the third line of Eq. (16.10 . which depends on the Laplacian 
of the potentials. This term represents gravitational lensing and it agrees with an existing 
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expression for the convergence in the context of modified gravity theories [24]. Consider 
next the second term in the third hne, which contains two time derivatives. The fact that 
this term is integrated with the standard lensing weight function (x — Xo) {Xs ^ x)/ iXs — Xo) 
signals that this term is also a lensing term. However, it is smaller than the leading lensing 
term by ~ (aH/k)'^, where k denotes the comoving wave number of a perturbation, and H 
denotes the Hubble rate. Therefore, it is subdominant for subhorizon perturbations. 

There is some ambiguity in the above classification of the various terms in Eq. (16.11) . 
and some of the terms may be classified differently. For example, we may use integration by 
parts to replace some of the boundary terms with terms containing integrals over derivatives 
of the potentials. The ambiguity in the classification originates from the fact that only the 
entire luminosity distance is observable, and the different individual terms do not correspond 
to a gauge invariant expression. 

As mentioned in the introduction, in many cases Eq. (16.11) can be simplified by spe- 
cializing to subhorizon perturbations. Under this approximation, terms that have a time 
derivative of the potentials are smaller by ~ aH/ k with respect to the term containing spa- 
tial derivatives. Furthermore, under this approximation Eq. (II. 3p implies that the potential 
ip is smaller by ~ aH/ k with respect to the velocity terms. While the potential is uncon- 
strained by relation (II. Sp . it is reasonable to assume that its magnitude is not much larger 
then the magnitude of ip. For this reason we neglect terms that depend on (p and ip (with 
no derivatives). The above considerations gives the approximate expression 

Dr\z, n) ^ (;^, - ;^„)(1 + ;.)|l + V, ■ n - , " "^^^ 

iXs - Xo)T-Cs 

^ Jxo Xs Xo "* 

that was discussed in the introduction. 
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